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Abstract 

Extending the BCS theory towards the strong electron-phonon interaction (EPI), a charged Bose 
liquid of small bipolarons has been predicted by us with a further prediction that the highest su- 
perconducting critical temperature is found in the crossover region of the EPI strength from the 
BCS-like to bipolaronic superconductivity. Later on we have shown that the unscreened (infinite- 
range) Frohlich EPI combined with the strong Coulomb repulsion create superlight small bipo- 
larons, which are several orders of magnitude lighter than small bipolarons in the Holstein-Hubbard 
model (HHM) with a zero-range EPI. The analytical and numerical studies of this Coulomb- Frohlich 
model (CFM) provide the following recipes for room-temperature superconductivity: (a) The par- 
ent compound should be an ionic insulator with light ions to form high-frequency optical phonons, 
(b) The structure should be quasi two-dimensional to ensure poor screening of high-frequency 
phonons polarized perpendicular to the conducting planes, (c) A triangular lattice is required in 
combination with strong, on-site Coulomb repulsion to form the small superlight bipolaron, (d) 
Moderate carrier densities are required to keep the system of small bipolarons close to the Bose- 
Einstein condensation regime. Clearly most of these conditions are already met in the cuprates. 

PACS numbers: 71.38.-k, 74.40. +k, 72.15.Jf, 74.72.-h, 74.25.Fy 



I. INTRODUCTION 

The discovery of high-temperature superconductivity in cuprates [1] has widened signif- 
icantly our horizons of the theoretical understanding of the physical phenomenon. A great 
number of observations point to the possibility that the cuprate superconductors may not be 
conventional Bardeen-Cooper-Schrieffer (BCS) superconductors [2J, but rather derive from 
the Bose-Einstein condensation (BEC) of real-space small bipolarons jjj, 4, 5L Importantly 
a first proposal for high temperature superconductivity by Ogg Jr in 1946 6], already in- 
volved real-space pairing of individual electrons into bosonic molecules with zero total spin. 
This idea was further developed as a natural explanation of conventional superconductivity 
by Schafroth [z| and Butler and Blatt Unfortunately the Ogg-Schafroth picture was 
practically forgotten because it neither accounted quantitatively for the critical behavior of 
conventional superconductors, nor did it explain the microscopic nature of attractive forces 
which could overcome the Coulomb repulsion between two electrons constituting a pair. 
On the contrary highly successful for low-T c metals and alloys the BCS theory, where two 
electrons were indeed correlated, but at a very large distance of about 10 3 times of the 
average inter-electron spacing, led many researchers to believe that any superconductor is a 
"BCS-like". 

However it has been found — unexpectedly for many researchers — that the BCS theory 
and its extension 9J towards the intermediate coupling regime, A < 1, break down already 
at A > 1 10j. It happens since the Migdal " noncrossing" approximation 11[ of the theory 
is not applied at A > 1. In fact, the small parameter of the theory, Xuq/Ef, becomes large 
at A > 1 because the bandwidth is narrowed and the Fermi energy, Ep is renormalised 
down exponentially due to the small polaron formation [h], [3] (here u>o is the characteristic 
phonon frequency, and we take h = c = ks = 1)- Extending the BCS theory towards the 
strong interaction between electrons and ion vibrations, A>1, a charged Bose gas of tightly 
bound small bipolarons was predicted [13;] instead of Cooper pairs, with a further prediction 
that the highest superconducting transition temperature is attained in the crossover region 
of EPI strength, A ~ 1, between the BCS and bipolaronic superconductivity jioj ] . 

For a very strong EPI polarons become self-trapped on a single lattice site and bipolarons 
are on-site singlets. A finite on-site bipolaron mass appears only in the second order of po- 



laron hopping, [13], so that on-site bipolarons might be very heavy in HHM, where EPI is 
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short-ranged. Actually HHM led some authors to the conclusion that the formation of itin- 
erant small polarons and bipolarons in real materials is unlikely [14] , and high-temperature 
bipolaronic superconductivity is impossible 15|. Nevertheless treating the onsite repulsion 
(Hubbard U) and the short-range EPI on an equal footing led several authors to the opposite 
conclusion with respect to bipolaron mobility even in HHM, which is generally unfavorable 
for coherent tunnelling. Aubry 16| found along with the onsite bipolaron (SO) also an 
anisotropic pair of polarons lying on two neighboring sites (i.e. the intersite bipolaron, 
SI) with classical phonons in the extreme adiabatic limit. Such bipolarons were originally 
hypothesized in [171 . Il8| to explain the anomalous nuclear magnetic relaxation (NMR) in 
cuprate superconductors. The intersite bipolaron could take a form of a " quadrisinglet" 
(QS) in 2D HHM, where the electron density at the central site is 1 and "1/4" on the four 
nearest neigbouring sites. In a certain region of U, where QS is the ground state, the double- 
well potential barrier which usually pins polarons and bipolarons to the lattice depresses to 
almost zero, so that adiabatic lattice bipolarons can be rather mobile. 

Mobile 5*1 bipolarons were found in ID HHM using variational methods also in the 
non and near-adiabatic regimes with dynamical quantum phonons Ji], [2^. The intersite 
bipolaron with a relatively small effective mass is stable in a wide region of the parameters of 



HHM due to both exchange and nonadiabaticity effects 20(] . Near the strong coupling limit 
the mobile 51 bipolaron has an effective mass of the order of a single Holstein polaron mass, 
so that one should not rule out the possibility of a superconducting state of 5*1 bipolarons 
with s or rf-wave symmetry in HHM [19]. More recent diagrammatic Monte Carlo study 
2l| found 5*1 bipolarons for large U at intermediate and large EPI and established the 
phase diagram of 2D HHM, comprising unbound polarons, SO and SI domains. Ref. 21] 
emphasised that the transition to the bound state and the properties of the bipolaron in 
HHM are very different from bound states in the attractive (negative U) Hubbard model 



without EPI 



22]. 



In any case the Holstein model is an extreme polaron model, with typically highest possi 



ble values of the (bi)polaron mass in the strong coupling regime [23|, [24j, [25|, 126J. Many doped 
ionic lattices, including cuprates, are characterized by poor screening of high-frequency op- 
tical phonons and they are more appropriately described by the finite-range Frohlich EPI. 
The unscreened Frohlich EPI provides relatively light lattice polarons and combined with 
the Coulomb repulsion also "superlight" but yet small (intersite) bipolarons. In contrast 
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with the crawler motion of on-site bipolarons, the intersite-bipolaron tunnelling is a crab- 
like, so that the effective mass scales linearly with the polaron mass. Such bipolarons are 
several orders of magnitude lighter than small bipolarons in HHM [23f . Here I review a few 



analytical 



23 



23| and more recent Quantum Monte-Carlo (QMC) [28|, |29|] studies of CFM 
which have found superlight bipolarons in a wide parameter range with achievable phonon 
frequencies and couplings. They could have a superconducting transition in excess of room 
temperature. 

II. COULOMB-FROHLICH MODEL 

Any realistic theory of doped narrow-band ionic insulators should include both the finite- 
range Coulomb repulsion and the strong finite-range EPI. From a theoretical standpoint, 
the inclusion of the finite-range Coulomb repulsion is critical in ensuring that the carriers 
would not form clusters. The Coulomb repulsion, V c , makes the clusters unstable and lattice 
bipolarons more mobile. 

To illustrate the point let us consider a generic multi-polaron " Coulomb- Frohlich" model 
(CFM) on a lattice, which explicitly includes the finite-range Coulomb repulsion, V c , and the 
strong long-range EPI 23j, |27|]. The implicitly present (infinite) Hubbard U prohibits double 
occupancy and removes the need to distinguish the fermionic spin, if we are interested in the 
charge rather than spin excitations. Introducing spinless fermion operators c n and phonon 



operators d m , the Hamiltonian of CFM is written in the real-space representation as 27] 



H = T (n - n')c£cv + V c (n - n'^cj^ + (1) 

n^n' n^n' 

UJ ^2 ~ n )( e ' e m-n)4 C n(dJn + d m ) + 

m ^ ' 

where T(n) is the bare hopping integral in a rigid lattice. In general, this many-body model 
is of considerable complexity. However, if we are interested in the non or near adiabatic limit 
and the strong EPI, the kinetic energy is a perturbation. Then the model can be grossly 
simplified using the L aIlg -F, rS ov clonic, —nat, o» Q in the W aI1Mer _ tati o„ 
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for electrons and phonons, 



S = 22 9( m ~ n )( e " G m-n)c ] n c n {d ] m - d m ). 

Here we consider a particular lattice structure, where intersite lattice bipolarons tunnel 
already in the first order in T(n). That allows us to average the transformed Hamiltonian, 
H = exp(S)H exp(— S) over phonons to obtain 



H — Hq + H pert , 



where 



H = -E p ^ C n C n + ^ ^2 V ( U ~ n ') C n C nCn/C n / + UJ ^ ( d ln d m + ^ J 
n n^n' m ^ ' 



and 



Hpert = 2J t(n - n')c^Cn/. 
n^n' 

is a perturbation. i? p is the familiar polaron level shift, 

E p = uj ^2g 2 (m - n)(e ■ e m _ n ) 2 , 

mu 

which is independent of n. The polaron-polaron interaction is 

v(n - ri) = V c (n - ri) - V ph (n - n), 

where 

V ph (n - n) = 2^0*^2 g(m - n)g(m - n')(e • e m _ n )(e ■ e m _ n /). 

m 

The transformed hopping integral is t(n — n') = T(n — n') exp[— g 2 (n — n')] with 
g 2 (n-n) = #(m - n)(e • e m _ n ) x 

[g(m - n)(e • e m _ n ) - g(m - n')(e • e m _ n /)] 



(2) 



(3) 



(4) 



(5) 



(6) 



at low temperatures. The mass renormalization exponent can be expressed via E p and V p h 

as 



g (n- n 



1 

u 



1 . , ' 



(7) 



The Hamiltonian H, Eq.([2]), in zero order with respect to the hopping describes localised 
polarons and independent phonons, which are vibrations of ions relative to new equilibrium 
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positions depending on the polaron occupation numbers. Importantly the phonon frequen- 
cies remain unchanged in this limit at any polaron density, n. At finite A and n there is a 
softening of phonons 5ujq of the order of uj Q n/X 2 [31). Interestingly the optical phonon can 
be mixed with a low-frequ ency polaronic plasmon forming a new excitation, "plasphon", 



which was proposed in 3l|, ]32[] as an explanation of the anomalous phonon mode splitting 



observed in cuprates 33j. The middle of the electron band is shifted down by the polaron 
level-shift E p due to the potential well created by lattice deformation. 

When V p h exceeds V c the full interaction becomes negative and polarons form pairs. The 
real space representation allows us to elaborate more physics behind the lattice sums in V p h 



271 ] . When a carrier (electron or hole) acts on an ion with a force f, it displaces the ion by 



some vector x = f/k. Here k is the ion's force constant. The total energy of the carrier- 
ion pair is — f 2 /(2k). This is precisely the summand in Eq.([3]) expressed via dimensionless 
coupling constants. Now consider two carriers interacting with the same ion. The ion 
displacement is x = (fi + f 2 )/A: and the energy is —ff/(2k) — f|/(2/c) — (fi • f2)/k. Here 
the last term should be interpreted as an ion-mediated interaction between the two carriers. 
It depends on the scalar product of fi and f2 and consequently on the relative positions of 
the carriers with respect to the ion. If the ion is an isotropic harmonic oscillator, then the 
following simple rule applies. If the angle <fi between fi and f2 is less than 7r/2 the polaron- 
polaron interaction will be attractive, if otherwise it will be repulsive. In general, some ions 
will generate attraction, and some repulsion between polarons. 

The overall sign and magnitude of the interaction is given by the lattice sum in Eq.Q. 
One should note that according to Eq.(]7]) an attractive EPI reduces the polaron mass (and 
consequently the bipolaron mass), while repulsive EPI enhances the mass. Thus, the long- 
range EPI serves a double purpose. Firstly, it generates an additional inter-polaron at- 
traction because the distant ions have small angle <fi. This additional attraction helps to 
overcome the direct Coulomb repulsion between polarons. And secondly, the Frohlich EPI 



makes lattice (bi)polarons lighter. Here, following [23|, ]27|, ]28|, ]29] , we consider a few examples 
of intersite superlight bipolarons. 
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FIG. 1: Apex bipolaron tunnelling in perovskites (after 



m) 



III. APEX BIPOLARONS 

High-T c oxides are doped charged-transfer ionic insulators with narrow electron bands. 
Therefore, the interaction between holes can be analyzed using computer simulation tech- 
niques based on a minimization of the ground state energy of an ionic insulator with two 
holes, the lattice deformations and the Coulomb repulsion fully taken into account, but ne- 
glecting the kinetic energy terms. Using these techniques net inter-site interactions of the 
in-plane oxygen hole with the apex hole, Fig.l, and of two in-plane oxygen holes, Fig.2, 
were found to be attractive in La 2 CuO^ [34J] with the binding energies A = 119meV and 
A = 60meV, respectively. All other interactions were found to be repulsive. 

Both apex and in-plane bipolarons can tunnel from one unit cell to another via the 
smg/e-polaron tunnelling from one apex oxygen to its apex neighbor in case of the apex 



bipolaron 23[], Fig.l, or via the next-neighbor hopping in case of the in-plane bipolaron 271 ] . 



Fig.2. The Bloch bands of these bipolarons are obtained using the canonical transformation, 
described above, projecting the transformed Hamiltonian, Eq.([2]), onto a reduced Hilbert 
space containing only empty or doubly occupied elementary cells [4|. The wave function of 

the apex bipolaron localized, say in cell m is written as 

4 

|m) = ^Ac!4 pf JO>, (8) 
1=1 

where i denotes the p x>y orbitals and spins of the four plane oxygen ions in the cell, Fig.l, 
and cl pex is the creation operator for the hole in one of the three apex oxygen orbitals with 
the spin, which is the same or opposite to the spin of the in-plane hole depending on the 
total spin of the bipolaron. The probability amplitudes Aj are normalized by the condition 



FIG. 2: Four degenerate in-plane bipolaron configurations A, B, C, and D . Some single-polaron 
hoppings are indicated by arrows (Reproduced from A. S. Alexandrov and P. E. Kornilovitch, J. 
Phys.: Condens. Matter 14, 5337 (2002), (c) IOP Publishing Limited, 2007. [271). 

\Ai\ = 1/2, if four plane orbitals p x i,Py2,Px3 an d p y 4 are involved, or by \Ai\ = l/v2 if only 
two of them are relevant. Then a matrix element of the Hamiltonian Eq.Q describing the 
bipolaron tunnelling to the nearest neighbor cell m + a is found as 

t b = (m\H\m + a) = \A^T^ x e^\ (9) 

where T^ x e~ g2 is a single polaron hopping integral between two apex ions. The inter-site 
bipolaron tunnelling appears already in the first order with respect to the single-hole transfer 
T^ x , and the bipolaron energy spectrum consists of two subbands ^'^(K), formed by the 
overlap of p x and p y apex oxygen orbitals, respectively (here we take the lattice constant 
a = 1): 

E X (K) = tcos(K x ) -t'cos(K y ), (10) 
E y (K) = -t'cos(K x ) +tcos(K y ). 

They transform into one another under ir/2 rotation. If t,t' > 0, "x" bipolaron band has 
its minima at K =(±7r, 0) and y-band at K = (0, ±7r). In these equations t is the renormal- 
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ized hopping integral between p orbitals of the same symmetry elongated in the direction of 
the hopping (ppa) and if is the renormalized hopping integral in the perpendicular direction 
(ppn). Their ratio t/t' = T^, ex /T"^f x = 4 as follows from the tables of hopping integrals in 
solids. Two different bands are not mixed because Tff^, = for the nearest neighbors. A 

Px,Py U 

random potential does not mix them either, if it varies smoothly on the lattice scale. Hence, 
we can distinguish 'x' and ! y' bipolarons with a lighter effective mass in x or y direction, 
respectively. The apex z bipolaron, if formed, is ca. four times less mobile than x and y 
bipolarons. The bipolaron bandwidth is of the same order as the polaron one, which is a 
specific feature of the inter-site bipolaron. For a large part of the Brillouin zone near (0, it) 
for 'x' and (ir, 0) for ! y' bipolarons, one can adopt the effective mass approximation 

K 2 K 2 

E x 'y(K) = — ^- + — y — (ii) 

with K x ^ y taken relative to the band bottom positions and m** = 1/t, m** = 4m**. 

X and y bipolarons bose-condense at the boundaries of the center-of-mass Brillouin zone 
with K = (±7r, 0) and K = (0, ±7r), respectively, which explains the d-wave symmetry and 
the checkerboard modulations of the order parameter in cuprates 35 1. 



IV. IN-PLANE BIPOLARONS 

Now let us consider in-plane bipolarons in a two-dimensional lattice of ideal octahedra 
that can be regarded as a simplified model of the copper-oxygen perovskite layer, Fig. 3 j^ . 
The lattice period is a = 1 and the distance between the apical sites and the central plane 
is h — a/2 = 0.5. For mathematical transparency we assume that all in-plane atoms, both 
copper and oxygen, are static but apex oxygens are independent three-dimensional isotropic 
harmonic oscillators. 

Due to poor screening, the hole-apex interaction is purely coulombic, 

g a (m - n) = -, 

m — n\ z 

where a = x,y, z. To account for the fact that c axis-polarized phonons couple to the holes 



stronger than others due to a poor screening 
direct hole-hole repulsion is 



23 



24 



271 ]. we choose /y/2. The 



V, 



FIG. 3: Simplified model of the copper-oxygen perovskite layer. Reproduced from A. S. Alexandrov 
and P. E. Kornilovitch, J. Phys.: Condens. Matter 14, 5337 (2002), (c) IOP Publishing Limited, 

2007. bal 



so that the repulsion between two holes in the nearest neighbor (NN) configuration is V c . 
We also include the bare NN hopping Tjvjv, the next nearest neighbor (NNN) hopping 
across copper T NNN and the NNN hopping between the pyramids T' NNN . 

The polaron shift is given by the lattice sum Eq.(3), which after summation over polar- 
izations yields 

E p = 2k 2 x u ( , 1 ,4 + | h " , 6 ) = 31.15^o, (12) 

"I m-nr m-nV 

m \ i i i 

where the factor 2 accounts for two layers of apical sites. For reference, the Cartesian 
coordinates are n = (n x + l/2, n y + 1/2, 0), m = (m x , m y , h), and integers. 
The polaron-polaron attraction is 

* 9 s h 2 + (m — n') ■ (m — n) . „ s 

V ph n-n' )=4oj4 V r^- 13 

^— ' m - n' 3 m - n 3 

m 1 11 1 

Performing the lattice summations for the NN, NNN, and NNN' configurations one finds 
Vph = 1.23 E p , 0.80 E p , and 0.82 E p , respectively. As a result, we obtain a net inter-polaron 
interaction as vnn = V c — 1.23 E p , vnnn = ^% — 0.80 E p , v' NNN = % — 0.82 _E P , and the 
mass renormalization exponents as g 2 NN = 0.38(E p /lu), g 2 NNN = 0.60(E p /u) and (q'nnn) 2 = 
0.59{E p /u). 

Let us now discuss different regimes of the model. At V c > 1.23 E p , no bipolarons are 
formed and the system is a polaronic Fermi liquid. Polarons tunnel in the square lattice with 
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t = T NN exp(-0.38E p /uj) and t' = T NNN exp(-0.60£ p /cd) for NN and NNN hoppings, 
respectively. Since g 2 NNN ~ (q'nnn) 2 one can neglect the difference between NNN hoppings 
within and between the octahedra. A single polaron spectrum is therefore 

£i(k) = -E p -2t'[ cos k x + cos ky] — 4i cos(/c.j;/2) cos(fcj,/2). (14) 

The polaron mass is m* = l/(t + 2t'). Since in general t > t', the mass is mostly determined 
by the NN hopping amplitude t. 

If V c < 1.23 E p then intersite NN bipolarons form. The bipolarons tunnel in the plane 
via four resonating (degenerate) configurations A, B, C, and D, as shown in Fig. 2. In the 
first order of the renormalized hopping integral, one should retain only these lowest energy 
configurations and discard all the processes that involve configurations with higher energies. 
The result of such a projection is the bipolaron Hamiltonian, 

H b = (V c - 3.23 E p ) Y^AlAi + BlBx + CfCi + A f A] (15) 

l 

-H Y^AiBi + B l°i + C i A + D\A X + H.c.) 
1 

n 

where 1 numbers octahedra rather than individual sites, x = (1,0), and y = (0, 1). A Fourier 
transformation and diagonalization of a 4 x 4 matrix yields the bipolaron spectrum: 

E 2 (K) = V C - 3.23E P ± 2t'[cos{K x /2) ± cos{K y /2)). (16) 

There are four bipolaronic subbands combined in the band of the width 8t'. The effective 
mass of the lowest band is m** = 2/t'. The bipolaron binding energy is A m 1.23E p — V c . 
Inter-site bipolarons already move in the first order of the single polaron hopping. This 
remarkable property is entirely due to the strong on-site repulsion and long-range electron- 
phonon interactions that leads to a non-trivial connectivity of the lattice. This fact combines 
with a weak renormalization of t' yielding a superlight bipolaron with the mass m** oc 
exp(0.60 Ep/u). We recall that in the Holstein model m** oc exp(2E p /uj) [13]. Thus the 
mass of the Frohlich bipolaron in the perovskite layer scales approximately as a cubic root 
of that of the Holstein bipolaron. 

At even stronger EPI, V c < 1A6E P , NNN bipolarons become stable. More importantly, 
holes can now form 3- and 4-particle clusters. The dominance of the potential energy over 
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kinetic in the transformed Hamiltonian enables us to readily investigate these many-polaron 
cases. Three holes placed within one oxygen square have four degenerate states with the 
energy 2(V C - 1.23.E7p) + V c /V2 - 0.80E P . The first-order polaron hopping processes mix 
the states resulting in a ground state linear combination with the energy E 3 = 2.71V C — 
3.26E p — V 'At 2 + t' 2 . It is essential that between the squares such triads could move only 
in higher orders of polaron hopping. In the first order, they are immobile. A cluster of 
four holes has only one state within a square of oxygen atoms. Its energy is £4 = A(V C — 
1.23£p) + 2(V C /V2 - 0.80£p) = 5A1V C - 6.52E P . This cluster, as well as all bigger ones, are 
also immobile in the first order of polaron hopping. We would like to stress that at distances 
much larger than the lattice constant the polaron-polaron interaction is always repulsive, 
and the formation of infinite clusters, stripes or strings is prohibited. We conclude that at 
V c < 1.16Ep the system quickly becomes a charge segregated insulator. 

The fact that within the window, 1.16E P < V c < 1.23E P , there are no three or more 
polaron bound states, indicates that bipolarons repel each other. The system is effectively 
a charged Bose-gas, which is a superconductor {f3,Q]. This superconducting state requires a 
rather fine balance between electronic and ionic interactions in cuprates. 



V. ALL-COUPLING LATTICE BIPOLARONS 



The multi-polaron CFM model discussed above is analytically solvable in the strong- 
coupling nonadibatic (ujq > T(a)) limit using the Lang-Firsov transformation of the Hamil- 



tonian, Eq.(l), and projecting it on the inter-site pair Hilbert space 23|, |27|. To extend 



the theory for the whole parameter space an advanced continuous time QMC technique 



28 



29|. Using CTQMC refs. 



28 



29] 



(CTQMC) has been recently developed for bipolarons 
simulated the CFM Hamiltonian on a staggered triangular ladder (ID), triangular (2D) and 
strongly anisotropic hexagonal (3D) lattices including triplet pairing [29(. On such lattices, 
bipolarons are found to move with a crab like motion (Fig. 1), which is distinct from the 

n 

crawler motion found on cubic lattices [13|]. Such bipolarons are small but very light for a 
wide range of electron-phonon couplings and phonon frequencies. EPI has been modeled 
using the force function in the site-representation as 
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FIG. 4: Polaron to bipolaron mass ratio for a range of uo = uJo/T(a) and A on the staggered ladder. 
Mobile small bipolarons are seen even in the adiabatic regime Q = 0.5 for couplings A up to 2.5 
(Reproduced from J. P. Hague et al., Phys. Rev. Lett 98, 037002 (2007), (c) American Physical 



Society, 2007. 



28f ). 



H e _ ph - - ^ /m(n) cJ^CW^ . (17) 
nmcr 

Each vibrating ion has one phonon degree of freedom £ m associated with a single atom. 
The sites are numbered by the indices n or m for electrons and ions respectively. Operators 
c nCT creates an electron on site n with spin a. Coulomb repulsion V(n — n') has been screened 
up to the first nearest neighbors, with on site repulsion U and nearest-neighbor repulsion V c . 
In contrast, the Frohlich interaction is assumed to be long-range, due to unscreened inter- 
action with c-axis high-frequency phonons |23|] . The form of the interaction with c-axis po- 
larized phonons has been specified via the force function 24], / m (n) = k [(m — n) 2 + 1] 3 ^ 2 , 
where k is a constant. The dimensionless electron-phonon coupling constant A is defined as 
^ = E m /m(^/^ w2 ^( fl ) which is the ratio of the polaron binding energy to the kinetic 
energy of the free electron zT(a), and the lattice constant is taken as a = 1. 

In the limit of high phonon frequency uj ^> T(a) and large on-site Coulomb repulsion 
(Hubbard U), the model is reduced to an extended Hubbard model with intersite attraction 
and suppressed double-occupancy by applying the Lang-Firsov canonical transformation 
(section 2). Then the Hamiltonian can be projected onto the subspace of nearest neighbor 
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FIG. 5: Bipolaron radius (in units of a) for a range of u> and A on the staggered ladder (Reproduced 
from J. P. Hague et al., Phys. Rev. Lett 98, 037002 (2007), (c) American Physical Society, 2007. 



2J" 



intersite crab bipolarons (sections 4). In contrast with the crawler bipolaron, the crab 
bipolaron's mass scales linearly with the polaron mass (m** = Am* on the staggered chain 
and m** = 6m* on the triangular lattice). 

Extending the CTQMC algorithm to systems of two particles with strong EPI and 
Coulomb repulsion solved the bipolaron problem on a staggered ladder, triangular and 
anisotropic hexagonal lattices from weak to strong coupling in a realistic parameter range 
where usual strong and weak-coupling limiting approximations fail. Importantly small but 



light bipolarons have been founc 



phonon frequency, uj < T(a) (28 



for more realistic intermediate values of EPI, A < 1 and 



29]. 



Figure 4 shows the ratio of the polaron to bipolaron masses on the staggered ladder 
as a function of effective coupling and phonon frequency for V c = 0. The bipolaron to 
polaron mass ratio is about 2 in the weak coupling regime (A <C 1) as it should be for a 



large bipolaron 36]. In the strong-coupling, large phonon frequency limit the mass ratio 
approaches 4, in agreement with strong-coupling arguments given above. In a wide region of 
parameter space, we find a bipolaron/polaron mass ratio of between 2 and 4 and a bipolaron 
radius similar to the lattice spacing, see Figs. 4 and 5. Thus the bipolaron is small and light 
at the same time. Taking into account additional intersite Coulomb repulsion V c does not 
change this conclusion. The bipolaron is stable for V c < 4T(a). As V c increases the bipolaron 
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mass decreases but the radius remains small, at about 2 lattice spacings. Importantly, the 
absolute value of the small bipolaron mass is only about 4 times of the bare electron mass 
m , for A = uj/T{a) = 1 (see Fig. 4). 

Simulations of the bipolaron on an infinite triangular lattice including exchanges and 
large on-site Hubbard repulsion U = 20T(a) also lead to the bipolaron mass of about 
Smoxy an d the bipolaron radius Rb p ~ 2a for a moderate coupling A = 0.5 and a large 
phonon frequency to = T(a) (for the triangular lattice, mo X y — l/3a 2 T(a)). Finally, the 
bipolaron in a hexagonal lattice with out-of-plane hopping X" = T(a)/3 has also a light 
in-plane inverse mass, m** « 4.5mo xy but a small size, Rb p ~ 2.6a for experimentally 
achievable values of the phonon frequency u = T(a) = 200meV and EPI, A = 0.36. Out- 
of-plane m** « 70moz is Holstein like, where m 0z = l/2d 2 T', (d is the inter-plane spacing). 
When bipolarons are small and pairs do not overlap, the pairs can form a BEC at Tbec — 

If we choose realistic values for the lattice constants of 0.4 
nm in the plane and 0.8 nm out of the plane, and allow the density of bosons to be ne=0.12 
per lattice site, which easily avoids overlap of pairs, then Tbec ~ 300K. 



VI. SUMMARY 



For a very strong electron-phonon coupling in the Holstein model with the zero-range EPI, 
polarons become self-trapped on a single lattice site and bipolarons are on-site singlets. The 
on-site bipolaron mass appears only in the second order of polaron hopping [13] , so that 
on-site bipolarons are very heavy. This estimate led some authors to the conclusion that 
high-temperature bipolaronic superconductivity is impossible . 

However we have found that small but relatively light bipolarons could exist within the 
realistic range of the finite-range EPI with high-frequency optical phonons. The effect 
appears since the finite-range Frohlich interaction combined with the long-range Coulomb 
repulsion provides an effective interaction with a deep attraction minimum for two holes on 
the neighbouring sites, and repulsive for other hole configurations. Bipolarons which are 
both light and small give rise to Ogg-Schafroth's bose-condensed state of charged bosons 
at high-temperatures, since the Bose-Einstein condensate has transition temperature that 



is inversely proportional to mass. Our conclusion is backed up by analytical 23j, [271] and 



CTQMC studies [28 



These studies let us believe that the following recipes is worth 
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investigating to look for room-temperature superconductivity [28[ : (a) The parent compound 
should be an ionic insulator with light ions to form high-frequency optical phonons, (b) 
The structure should be quasi two-dimensional to ensure poor screening of high-frequency 
phonons polarized perpendicular to the conducting planes, (c) A triangular lattice is required 
in combination with strong, on-site Coulomb repulsion to form the small superlight crab 
bipolaron, (d) Moderate carrier densities are required to keep the system of small bipolarons 
close to the dilute regime. Clearly most of these conditions are already met in the cuprate 
superconductors. 
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